SQUARING THE CIRCLE AND CUBING THE SPHERE: 
CIRCULAR AND SPHERICAL COPULAS 

MICHAEL D. PERLMAN AND JON A. WELLNER 

Abstract. Do there exist circular and spherical copulas in M''? That is, do there 
exist circularly symmetric distributions on the unit disk in and spherically sym- 
metric distributions on the unit ball in K'^, d > 3, whose one-dimensional marginal 
distributions are uniform? The answer is yes for d = 2 and 3, where the circular 
and spherical copulas are unique and can be determined explicitly, but no for d > 4. 
A one-parameter family of elliptical bivariate copulas is obtained from the unique 
circular copula in by oblique coordinate transformations. Copulas obtained by 
a non-linear transformation of a uniform distribution on the unit ball in M'^ are also 
described, and determined explicitly for d — 2. 

KEY WORDS AND PHRASES: Bivariate distribution, multivariate distribution, 
unit disk, unit ball, circular symmetry, spherical symmetry, circular copula, spher- 
ical copula, elliptical copula. 



1. Introduction 

Do there exist spherically symmetric distributions on the closed unit ball in M*^ 
that have uniform one-dimensional marginal distributions on [—1, 1]? A distribution 
on Bd with this property may be said to "square the circle" when d = 2 and to "cube 
the sphere" when d > 3. 

The cumulative distribution function (cdf ) of a multivariate distribution on the unit 
cube [0,1]*^ whose marginal distributions are uniform [0,1] is commonly called a copula; 
see Nelsen (2006) for an accessible introduction to this topic. However, although it 
is customary to confine attention to distributions on the unit cube, our interest is 
in spherically symmetric (= orthogonally invariant) distributions on B^ with uniform 
marginal distributions. Therefore we take "copula" to mean a multivariate cdf on the 
centered cube Cd := [—1, l]'^ with uniform[-l,l] marginals. 

For d = 2 (resp., d> 3), such a copula, if it exists, will be called a circular copula 
(resp., spherical copula) if it is the cdf of a circularly symmetric (resp., spherically 
symmetric) distribution on the unit disk B2 (resp., unit ball Bd)- 

It will be noted in Sections |2] and [3] that circular and spherical copulas are unique if 
they exist, but exist only for dimensions d = 2 and d = 3. The proof of non-existence 
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for (i > 4 is remarkably simple. Explicit expressions for these copulas are given in 
Sections |3] and |4] respectively. 

In Section |5j a new one-parameter family of bivariate copulas called elliptical cop- 
ulas is obtained from the unique circular copula in by oblique coordinate trans- 
formations. Finally, in Section |6| copulas obtained by a non-linear transformation of 
a uniform distribution on the unit ball in are described, and determined explicitly 
for d = 2. 

2. Uniqueness and Existence of Circular and Spherical Copulas 

Proposition 2.1. Circular and spherical copulas are unique if they exist. 

Proof. If a circular or spherical copula exists on Cd, it is the cdf of a random vector 
Z = {Zi, . . . , Zd) with a spherically symmetric distribution on Bd and with each 
Zi ~ uniform [— 1, 1]. The latter implies that Z has no atom at the origin, i.e., 
P[Z = 0] = 0, so we may consider the "polar coordinates" representation Z = R ■ U, 



where R= \\Z\\ < 1 and U = Z/\\Z\\. It is well known (e.g., Cambanis et al. (1981), 
Lemmas 1 and 2) that the random unit vector U = {Ui, . . . ,Ud) is independent 
of R and is uniformly distributed on the unit sphere dBd, which implies that each 
Uf ~ Beta(l/2, {d - l)/2). Since Zi = RUi, we have that 

log(Z^) = log(/2^) + log(t/n. (2.1) 

Because R and Ui are independent, it follows that the characteristic function of 
log(i?^) is the quotient of the characteristic functions of \og{Zf) and log([/f). Thus 
the distribution of log(i?^), and therefore that of R, is uniquely determined by the 
distributions of Zf and U^, which are already specified above. Thus the the joint 
distribution of {R, U) is uniquely determined, hence so is the distribution of Z, hence 
so its cdf = copula. □ 

The existence of spherical copulas is easy to determine in three or more dimensions: 

Proposition 2.2. Spherical copulas do not exist for d > A. For d = 3, the unique 
spherical copula is generated by the uniform distribution on the unit sphere dB^ : = 
{{xi,X2,xs) \ xl + xl + xl = 1}. 



Proof. Let Z be as in the proof of Proposition 2.1 Then 



^ = EiZf) = EiR')E{U^) < i. (2.2) 

since Zi ~ uniform[-l, 1], < i? < 1, and Uf ~ Beta(l/2, {d - l)/2). Thus d < 3, 
so a spherical copula cannot exist when d> 4. 



^This result is well-known (e.g., Feller (19711, pp. 31-33, who uses "random direction" to indicate 



the uniform distribution of J7 G dB^), and reappears frequently (e .g. Arellano- Valle (2001), Theorem 



3.1). The essence of the result goes back at least to Schoenberg (|1938|. 
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Furthermore, if a spherical copula is to exist for = 3, it follows from (2.2) that 
its generating random vector Z G -B3 must satisfy E{R^) = 1, hence R = 1 with 
probability one. This can occur only if Z is uniformly distributed on the unit sphere 
dB^. But it is well knowiij^that this distribution does indeed have uniform marginal 
distributions on [—1, 1], hence generates the unique spherical copula for ci = 3. □ 

3. The Bivariate Case: the Unique Circular Copula 
The following three questions constitute an engaging classroom exercise. 

Question 1. Let {X,Y) be a random vector uniformly distributed on the unit disk 
(= ball) B2 in M^. Find the marginal probability distributions of X and Y. 

Answer 1. One can easily show that X has the "semi-circular" probability density 
function (pdf) given by 

2 



/(x) = -Vl-x2, -1<X<1. (3.1) 
vr 

(See Figure [1]) By symmetry, Y has the same pdf as X. 

Question 2. Let (X, Y) be a random vector uniformly distributed on the unit circle 
dB2 in ]R2. Find the marginal probability distributions of X and Y. 
Answer 2. We can represent {X,Y) as (cosG, sinG) where G ~ uniform [0, 27r). It 
follows readily that X has pdf 

fix) = Jl ^ , -1< X < 1. (3.2) 



vr 



(See Figure [T]) By symmetry, Y has the same pdf as X. 




Figure 1. The densities (3.1) (lower, blue) and (3.2) (upper, purple) 



'This follows 
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In both cases, the joint distribution of {X, Y) is circularly symmetric, that is, 
invariant under all orthogonal transformations of M^. A comparison of the shapes of 
the pdfs in Figure [T] suggest a third question: 

Question 3. Does a circularly symmetric bivariate distribution with uniform[— 1, 1] 
marginals exist on i?2? If so, it determines a circular copula on C2, which is unique 
by Proposition 2.1 This also follows from uniqueness results for the Abel transform; 
see e.g. Bracewell (1986). 

Answer 3. Optimistically, let's seek an absolutely continuous solution. That is, we 
seek a bivariate pdf on B2 of the form 

such that the marginal pdf 

/\/l—x'^ 
f{x,y)dy, -1 < X < 1, 

is constant in x. Here g is a nonnegative function on (0, 1) that must satisfy 

2n [ rg{r^)dr = 1 (3.3) 



in order that ff^^ f{x,y)dxdy = 1 (transform to polar coordinates: {x,y) — )■ {r,9)). 
To determine a suitable g, first set h{t) = g{l — t), then let u = ^^f^ to obtain 

fix) = hil-x'-y')dy 

= [ h{{l - u^){l - x^))du 



2VI - X2 [ 

Jo 



h{{l ~ u^){l ~ x^))du. 



If we take h(t) = ct then clearly f{x) does not depend on x, and choosing c = l/27r 
satisfies (3.3). Thus the bivariate pdf 

/(x,i/) = — = ^ x' + y'<l, (3.4) 

ztta/ 1 — — 



determines a circularly symmetric bivariate distribution on B2 and yields the desired 
circular copula. 



Question 4. Having determined the unique circularly symmetric distribution (3.4) 
on B2 with uniform marginals, what is the corresponding cdf F{x,y), that is, what 
is the corresponding circular copula? 

4 



Circular and Spherical Copulas 



January 4, 2011 




Answer 4. The circular symmetry of {X, Y) implies that its distribution is invariant 

under sign changes, i.e., {X,Y) = (±X, ±y). By the following lemma, the cdf 
F{x,y) = P[X < X, Y < y] on C2 = [—1, 1]^ can be expressed in terms of Fo(x,?/), 
its truncation to the first quadrant: 

Fo(x, y) = P[0<X<x,0<Y<y] (3.5) 

for < x,y < 1, and also in terms of the complementary cdf F{x, y) = P[X > x, Y > 
y] for < x,y < 1. Because {X,Y) = (±X, ±y) and has uniform[— 1, 1] marginals, 

Fo(x, y) = P[0 < X < 1, < r < 1] - P[X > X, < F < 1] 
-P[0 <X <1,Y >y] + P[X >x,Y >y] 

F{x,y) 



4 V 4 y V 4 

= + 0<x,y<l. (3.6) 

Lemma 3.1. Let {X,Y) be a bivariate random vector on C2 with uniform[— 1,1] 

marginal distributions and sign-change invariance, i.e., (X,y) = (±X, ±Y). Then 
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for {x,y) e C2, 

F{x,y) 



X + y + 1 
4 

X + y + 1 



+ (^{xy) Fo{\x\, \y\) 



a{xy) 



\x\ + \y\-l 



n\AAy\) 



(3.7) 
(3.8) 



where a{w) = sign(-u;) if w and (t(0) = 0. 
Proof. To obtain ( |3.7 ), consider four cases: 



Case 1: < x,y < 1. Because (X, Y) is sign-change invariant and has uniform[— 1, 1] 
marginals, 

F{x, y) = P[0 < X < X, < Y < y] + P[0 < X < X, Y < 0] 

+p[x < 0, < r < + p[x < 0, r < 0] 

^ / „ X y 1 
= Fo{x,y) + - + '- + - 



X + y + 1 
4 



{xy)Fo{\x\, \y\). 



Case 2: —l<x<0<y<l. Similarly, 

F{x,y) = P[X <0,0<Y <y]-P[x <X <0,0<Y <y] 

+p[x < 0, r < 0] - p[x < X < 0, r < 0] 

= |-Fo(-a:,,) + --^ 



X + y + 1 



+ a{xy) Fo{\x\, \y\). 



F{x,y) 



Case 3: — 1 < y < < x < 1. Similarly, 

= P[0 < X < a;, r < 0] - P[0 < X < a;, ?/ < r < 0] 

+p[x < 0, r < 0] - p[x < 0, y < F < 0] 

X ^ , .1 (— y) 
= i-^o(x,-y) + --^ 

= ^ + a(x?/)Fo(F| 
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Case 4-' —^^x,y<0. Similarly, 



F{x, y ) = P[X < 0, F < 0] - P[a; < X < 0, r < 0] 

-P[X <0, y <Y <0] + P[x < X <0, y <Y <0] 

1 (-^) i-y) 

4 4 



+ Foi-x, -y) 



X + y + \ 



a{xy) Fq{\xI \y\). 



Finally, (3.8) follows from (3.7) by (3.6) 



□ 



Thus, to determine the circular copula F{x,y) for the pdf (3.4), it suffices to 
determine the complementary cdf F{x,y) for < x,y < 1 and apply (3.8). Because 
F{x,y) = when + > 1, we need only consider the case where + y"^ < 1. 



First approach: When < x,y < 1 and x^ + y'^ < 1, F{x,y) can be expressed as 
follows. By using Figure [3] we find that 




Figure 3. Region of integration, 2-dimensional case 
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F{x,y) 




1 

2^ 



arcsm 



ds. 



(3.9) 



However, we were unable to evaluate this integral directly. 



Second approach: Fortunately, we have found a solution in the molecular biology and 
optics literatures, where the problem of finding the area of the intersection of two 
spherical caps on the unit sphere dB^ has been addressed. The following general 



result is due to Tovchigrechko and Vakser (2001) and also appears in Oat and Sander 
( |2007| ). 

Lemma 3.2. Let Si and S2 be spherical caps on dB^. Let ri and r2 denote their 
angular radii and let d denote the angular distance between their centers (0 < d < tc). 
Assume that < ri,r2 < tt/2 and d < ri + r2, so that the intersection 5*1 fl 5*2 7^ 
and consists of a single "diangle"; (see Figures^ and^ ) Then Area^SiH S2) is given 

by 



A{ri,r2; d) 



27r 



2tx cos(ri) — 2n cos(r2) — 2 arccos 

cos((i) cos(ri) — cos(r2) 



cos((i) — cos(ri) cos(r2) 



smirij smlr2 1 



+ 2 cos(ri) arccos 
+ 2 cos(r2) arccos 



sin((i) sin(ri) 
cos((i) cos(r2) — cos(ri) 



sin((i) sin(r2) 

This result can be applied to obtain our desired circular copula as follows. 



(3.10) 



If (X, y, Z) is uniformly distributed on OB-^, then the event {X > x^Y > y} 
corresponds to the intersection of the two spherical caps {X > x} and {Y > y}, so 
P[X > X, Y > y] is given by the area A{x, y) of this intersection divided by the total 
area of dB^^ i.e., by An. (See Figures |4] and [sj) Also, the joint distribution of {X,Y) 
is circularly symmetric on the unit disk B2 and has uniform marginals, so must be 
the unique such bivariate distribution, namely the distribution with pdf (3.4). 
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Figure 4. Intersection of two spherical caps 




Figure 5. Intersection of two spherical caps, circular representation 
(modified from Tovchigrechko and Vakser (2001)) 
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Thus, for < x,y < 1 and + < 1, our desired complementary cdf is given by 

1 



F{x,y) 



An 
1 



A{x,y) 

y4(arccos(x), arccos(?/); 7r/2) 



An 

1 X y 1 
2~ 2~ 2~ Jn 



X 

H arccos 

2tx 



(3.11) 
(3.12) 



arccos 



xy 



x^ 



v/(l-x2)(l-y2) 

y 



—X 



2-n 



arccos 



1 — X — y 
4 



+ a(a;,y), 



(3.13) 



where for < x, ?/ < 1 and x^ + < 1, 



a{x,y) 



1 

2^ 



X arcsm 



arcsm 



x^ 



+ ?/ arcsin 



X 



xy 



Theorem 3.3. The unique circular copula on C2 is given by 

X + y + 1 



F{x,y) 



+ a{x,y), 



where a{x,y) is defined by (3.14) for x^ + < 1 and by 

a{x,y) = a[xy) 



(3.14) 



(3.15) 



(3.16) 



for x^ + 2/^ > 1. Note that (3.14) and (3.16) agree when x^ + = 1 and both are 
sign-change equivariant on C2: for all {x,y) G C2 and all e,6 = ±1, 

a{ex,6y) = e6 ■ a{x,y). (3-17) 



Proof. From (3.8) and (3.13), when x^ + < 1 we have 

X + y + 1 



F{x,y) 



X + y + 1 



+ a{xy)a{\x\, \y\) 
+ o:{x,y) 



(3.18) 
(3.19) 



by ([3l7l ). When x^ + > 1^ |^|) = so ( |3l5| ) again holds by (Q and 

(|3.16D. □ 

See Figure [6] for a plot of the resulting copula (on [—1, 1]^). 
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Figure 6. Copula (joint distribution function), Theorem 3.3 



4. The Trivariate Case: the Unique Spherical Copula 

Question 5. Having determined the unique spherically symmetric distribution on 
i?3 with uniform marginals, namely, the uniform distribution on the unit sphere dB^, 
what is the corresponding cdf F{x,y,z) on C3, i.e, the unique spherical copula? 



Answer 5. As in Section |3| let (A, F, Z) be uniformly distributed on dB^, so that 
F[x,y, z) = P[X < X, Y < y, Z < z]. Again we first determine the complementary 
cdf F{x, y, z) = P[X > X, Y > y, Z > z] ioi < x,y, z < 1 and + + < 1, 
the intersection of the first octant of C3 with the interior of B3. Here the event 
{A > X, Y > y, Z > z} corresponds to the intersection of the three spherical caps 
{A > x}, {Y > y}, and {Z > z} on dB^, so F{x,y,z) is the area A{x,y,z) of this 
intersection divided by the total area Atc of dB^. 

Recall that two approaches were proposed in Section |3] to obtain the area A{x,y) 
of the intersection of two circular caps {A > x} and {Y > y}. The first approach 



led to the integral (3.9) that we were unable to evaluate explicitly, so we adopted 



a second approach based on the geometric Lemma 3.2 of Tovchigrechko and Vakser 
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(2001). Audrey Tovchigrecliko has kindly suggested a method for extending Lemma 



3.2 to the case of three spherical caps in general position, which if carried out would 



yield an explicit expression for A(x,y,z). However, we have found that because the 
axes of our three caps are mutually orthogonal, the two approaches just mentioned 
for the bivariate case can be combined to obtain F{x,y, z) = ^A{x,y, z) directly for 
the trivariate case, as now described. 



We begin by extending (3.9) to obtain an integral expression for F(x,y,z) when 
< a;, y, 2 < 1 and x"^ + y"^ + z'^ < 1. We require the fact that 



< a, 6 < 1 and + 6^ = 1 implies arcsin(a) + arcsin(6) = n/2 
Lemma 4.1. If < x,y,z < 1 and x"^ + y"^ + z"^ < 1, then 

TT . f y 



F{x,y,z) 



1 

4:71 



— arcsm 



arcsm 



(4.1) 



ds. (4.2) 



Because (X, Y, Z) is exchangeable, (4.2) remains valid under any permutation ofx, y, z 
on the right-hand side. 



Proof. Since + + 
using Figure [7] that 



1 and {X,Y,Z) = {X,Y,-Z), it follows from (3.4) by 



P[X > x,Y > y, Z > z] 
= Ip[X >x,Y>y,X^ + Y^<l-z^] 



1 

Air 



VI - S2 - t2 



dt > ds 



1 

Air 
1 

Att 
1 

Air 



Vl-s2-22 



dt 



ds 



dv 



ds 



arcsm 



1 - s2 - ^2 

l-s2 



arcsm 



ds. 



Now apply (4.1) to obtain (4.2). 



□ 



As noted above, the integral in (4.2) appears difficult to evaluate explicitly, but 
the following indirect argument succeeds. Recall from (3.9) and (3.13) that when 
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Figure 7. Region of integration, 3-dimensional case, Lemma 4.1 



< x,y < 1 and + y'' < 1 



F{x,y) 



1 rV^y 

1 — X — y 



■ y 

arcsm , 

2 \VT^ 



ds 



+ a{x,y), 



where a{x,y) is given by (3.14). Because z < a/1 — x"^ — y"^ < a/1 — when < 
x,y,z < 1 and + + < 1, it follows that 



1 

2^ 



y/l-x^-y^ 



arcsm 



ds 



1 — x^ ~ y^ — z 



13 



(4.3) 
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Therefore from (4.2) and (4.3), ii < x,y,z < 1 and x + y + z < 1 then 



AttF^x, y, z) 



A/l-a;2-y2 r 



TT fx 

arcsm 



ds 



■ y 

arcsm , 

2 Vv^r^ 



ds 



IT 



1 — x^ — y^ — z 



— — x'^ — y'^ — z) + 2tt a{z, x) — a{\/l — x'^ — y'^, x) 



+ 2 (v^^ - x2 - t/2 _ 2) + 27r a{z, y) - a{y^ 1 - x"^ - y"^, y) 



(a/1 - - y2 _ ^) 



7r 



1 — — ?/2 — 2; 



+27r 



x) + q;(2;, — a{\Jl — x"^ — y"^, x) — ^(a/I — x^ — y^, y) 



(a/1 - X^ - ?/2 _ 2;) 



. , \ .IX 

+ X arcsm , + 2; arcsm 



X2 



y arcsm | — , | + z arcsm 



arcsm 



arcsm 



'1 — x^ — y2 



X arcsm 



+ arcsin 



1 — x^ — arcsin 



v/(l-a;2)(l-^2) 



v/(l-y2)(l-z2) 

X 



Xa/F— x2 — t/2 



v/(l -a;2)(x2 + |/2) 
2/ 



x2 — ?/2 arcsin 



y arcsm 



+ arcsin 



1 — x2 — ?/2 

7r 



^a/I - x2 - y2 



By (4.1), however. 
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so if we define h{x, y) by 

xy \ , ■ f a; a/1 - - y2 

h[x,y) = arcsm = + arcsm ' 



v/(l-x2)(l-?/2) / I v^r^^Va;2 + y2 



. , ?/Vl-x2-y2 
+ arcsm — ^^^^ (4.4) 



^/Y^^^/x^^^ 
ioT < x,y < 1 and + < 1, then 
ATTF{x,y,z) 

= — + — arcsm 



2 ^^(l_^2)(l_^2) 

+ X arcsin ( , ] + z arcsin ( , I — arcsin 

+ y arcsin I — -^^=^= ) + z arcsin I ^ I — arcsin 




\VT^J \y/{l-y^){l-z^) 



1 — x2 — y2 \ ^ / ^1 _ _ ^2 



X arcsm , — — y arcsm , , 

TT 

--Z + + a{x,y) 

■ f y \ .fx 

— X arcsm , — y arcsm — , 

z \ fx 

+ X arcsin I , + z arcsm , I — arcsin 

+ y arcsin I — , ] + z arcsin ( , ^ I — arcsin , , 

' KVi^J \v(i-t/2)(i-^2) 




1 — x2 — y2 \ _ / _ 3;2 _ ^2 



— X arcsm - — ^ — — y arcsm 



„,2 



where a{x,y) is given by (3.14). Now (4.1) gives 



- , y V . / Vl-2;2-y2\ 

X arcsm , + x arcsm — — , — ] = x { — 



X \ . / a/I -a;2 -y2 \ 
V arcsm — , + y arcsm , — \ = y 
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so the above simplifies to 



where 



AnF{x, y,z) = --{x + y + z) + h{x, y) + 27rA(a;, z), 



A(a;, 2;) = + a(a;, z) + 2;), 



(4.5) 
(4.6) 



a symmetric function of {x,y^z). By (4.1), however 
/i(x,?/) 



arcsm 



a/1 - x2 - ?/2 



v/(l-x2)(l-y2) 



+ arcsm 



xa/1 — a;2 — 



+ arcsin 

TT 



V^l -^/^V'x^ + yS 
7 + a + /3, 



and 



sin(a + j3) 
= sin a cos (3 + cos a sin (3 



X 



a/1 - a;2 _ y2 



a; 



+ 



l/V^l - X2 - y2 



a/1 — x2a/x2 + ?/2 y/l — y2^2;2 + 1/2 a/1 — a;2 ^a;2 + 7/2 y^l — 7/2 y^a;2 + 1/2 

^ a/1 -X2 - 7/2 

V(l-a;2)(l-7/2) 
= sin 7, 

so a + /3 = 7, hence /;.(x, y) = f identically in {x,y). Therefore we conclude that 

1 — X — 7/ — 2; A(x,7/,2;) 



F{x,y,z) 

for < a;,7/, z < 1 and + 7/^ + 2;^ < 1. 



(4.7) 



We now apply (4.7) to obtain the cdf F{x,y,z) for all {x,y,z) G C3. For this, 
extend the definition of A in (4.6) to all {x,y, z) e C3 by means of (3.14) and (3.16). 

Theorem 4.2. The unique spherical copula F{x,y,z) on C3 is given as follows: 
for + 7/^ + < 1, 



F{x,y,z) 



l+x+y+z _|_ A{x,y,z) 

8 ~'~ 2 ' 

l+x+y+z I A{x,y,z) 

8 ' 2 



+ a{xyz) 



l-\x\-\y\-\z\ A(|x|,|y|,|^|) 
8 ~'~ 2 

16 



if x^ + 7/^ + 2:^ < 1; 
if a;^ + 7/2 + 2^ > 1. 
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Proof. We use repeatedly the facts that (X, Y, Z) = (±X, ±Y, ±Z) and (X, Y, Z) 
has uniform[— 1, 1] marginals. First, by (3.15) and ( |3.17[ ), the complementary two- 
dimensional marginal cdf is given on all of C2 by 



— 1 — X — y 

F{x,y) = F{-x,-y) = ^a{x,y), {x,y) e C2. (4. 



Case la: < x,y, z, x + y + z < 1. By inclusion-exclusion, (4.8 ), and (4.7) 



F{x,y, z) 

= 1 - P[X >x\- P\Y >y]- P[Z > z] + P[X >x,Y >y]+ P[X > x, Z > z] 
+ P[Y >y, Z > z]- P[X > x,Y > y, Z > z] 

1 — X I — y 1 — z 1 — X — y , , 1 — X — z , , 
1 : : : 1 : 1- c([x, y) H h <y[x, zj 



4 

1-y-z 1-x-y-z A{x,y,z) 
+ ^^ + aiy,z) — 

1 + x + y + z A{x,y,z) 



4 



Case lb: < x,y, z, x"^ + y'^ + z"^ > 1. Here P[X > x, Y > y, Z > z] = 0, so 



F{x,y, z) 

= 1 - P[X >x]- P[Y >y]- P[Z > z] + P[X >x,Y>y]+ P[X > x, Z > z] 
+ P[Y >y, Z > z] 

1 — X 1 — y 1 — z 1 — X — y 1 — X — z 
1 : : : \ : h a{x, y) H h a(x, z) 



2 2 2 4 

I — y — z . , 
+ — ^ — + a{y,z) 

l + x + y + z A{x,y,z) 

^ + — ^ — + ^w) 

17 
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Case 2a: x < < y, z, + y"^ + < 1. By inclusion-exclusion, (4.8), and (4.7), 

F{x,y,z) 

= P[X <x]- P[X <x,Y >y]- P[X <x, Z > z] + P[X <x,Y>y,Z>z] 
^ ^ ^ P[-X <x,Y>y]- P[-X <x, Z> z] + P[-X <x,Y>y,Z>z] 

P[X > -X, Y >y]- P[X >-x,Z>z] + P[X > -x,Y > y, Z > z] 

x+1 1+x—y 1+x—z 

— a{—x, y) a(— X, z) 



2 

X + 1 



2 4 

^ l+x-y~ z _^ A(-x,y, z) 



4 



1 + x + y + z^ A{x, y, z) 



Case 2b: x < < y, z, x"^ + y'^ + z"^ > 1. Here P[X > -x, Y > y, Z > z] = 0, so 
F{x,y,z) 

= P[X <x]- P[X <x,Y >y]- P[X <x, Z> z] 
= - P[X >-x,Y>y]- P[X >-x,Z> z] 



x + 1 1 + X — y 



- a{-x,y) 



1 + x — z 



2 4 
1+x + y + z^ A{x, y, z) 



+ a{xyz) 



— «(— x, z) 

A - \y\ - \A , \yl \z\) 



Case 3a: y < < x, z, x^ + y"^ + z^ < 1. Similar to Case 2a. 

Case 3b: y < < x , z , x"^ + y"^ + z'^ > 1. Similar to Case 2b. 

Case 4-a: z < < x,y, x^ + y"^ + z^ < 1. Similar to Case 2a. 

Case 4b: z < < x,y, x"^ + y'^ + z"^ > 1. Similar to Case 2b. 

Case 5a: x,y < < z, x^ + y"^ + z^ < 1. By inclusion-exclusion, (4.8), and (4.7), 

F{x,y,z) = P[X <x,Y <y]- P[X <x,Y <y, Z > z] 

= P[X >-x,Y> -y] - P[X >-x,Y> -y, Z > z] 

1 + x + y 1 + x + y — z A(—x,—y,z) 
= -^— + a{-x,-y) 

1 + x + y + z A{x,y,z) 
8 ^ 2 ■ 
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Case 5b: x,y < < z, + y"^ + > 1. Here P[X >-x,Y> -y, Z > z]=0,so 
F{x,y,z) 
= P[X <x,Y <y] 
= P[X >-x,Y> -y] 
1 + X + y 



\x\ 



\y\ 



+ a{-x,-y) 
1 + x + y + z A{x,y,z) 

o ^ n ^ (^{xyZ) 



Case 6a: x , z < Q < y , x"^ + y"^ + z"^ < 1. Similar to Case 5a. 
Case 6b: x,z<0<y,x'^ + y'^ + z'^> 1. Similar to Case 5b. 
Case la: 2; < < x, + + 2;^ < 1. Similar to Case 5a. 
Case lb: y , z < Q < x, x"^ + y"^ + z"^ > 1. Similar to Case 5b. 

Case 8a: x,y, z < < z, x'^ + y'^ + z'^ < 1. 

F{x, y, z) = P[X < x,Y <y, Z < z] 

= P[X > -x,Y > -y, Z > - 
1 + x + y + z^ A{—x, —y, 



1 + x + y + z^ A{x, y, z) 



Case 8b: x,y,z < 0, x"^ + y"^ + z"^ > 1. 
F{,x,y,z) = 

1 + x + y + z A{x,y,z 



+ 



+ a{xyz) 



\x\ 



\y\ - 



\y\, \z\) 



□ 



5. A One-parameter Family of Elliptical Copulas 



Let (X,Y) ~ f{x,y) in (3.4), the unique circularly symmetric distribution on the 
unit disk B2 with uniform[— 1, 1] marginals. For any angle 7 G (— 7r/2, 7r/2), consider 
the transformed variables 



U = X, VI = Xsin7 + FC0S7. 



(5.1) 



By the circular symmetry of (X, Y), = y ~ uniform [—1, 1], so the random vector 
{U, V^) again generates a copula on the centered square C2. Denote the pdf and cdf 
of {U,Vy) by f^{u,v) and F^{u,v) respectively. Then {F^ \ 7 G (— 7r/2, 7r/2)} is a 
one-parameter family of elliptical copulas, so-called because the support of {U, Vy) is 
the ellipse 

:= {(m, v) \ + v"^ — 2nf sin7 < cos^7}. (5.2) 
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(Note that Eq = i?2-) From (5.1), the correlation coefficient of U and V.y is given 
simply by 

p(t/,K,) = sin7, (5.3) 
so 7 indicates the degree of linear dependence between U and V^. 

Proposition 5.1. The pdf of {U, Vy) is given by 

1 



-AeJu,v). 



(5.4) 



2n a/ cos^ 7 — — f ^ + 2uv sin 7 
Proof. The pdf can be obtained by a standard Jacobian computation. From ( |5.1 ), 

(5.5) 



V — usiwy 

x{u, V) = u, y{u, V) = —7777^, 



so 



dx{u,v) 



dx{u,v) 



cos(7) 



0, 



^fe--tan7, MM-^ 



du 



dv cos 7 ■ 



Thus the Jacobian of the transformation is J = 1/ cos 7, so from (3.4) we obtain 



fj{u,v) = f{x{u,v),y{u,v)) ■ J 

1 



l-K^X — V? — {v — usin7)^/ cos^ 7 

cos 7 1 

27r ^/(l - m2)cos2 7 - {y - u sin 7)2 cos 7 

1 , , 

=1^^(m, f ) 

27ry cos^ 7 — — + 2mw sin 7 



\B^{x{u,v),y{u,vy) 



cos 7 

lB2(x(M,f),t/(M,f)) 



□ 



Figure |8] shows the density f^^iu.v) with 7 = 7r/4. 

To describe the family of elliptical copulas F^, we extend the definitions (3.14) and 



(3.16) as follows. First, for {u,v) G -E^ define 

V — u sin 7 



1 
2^ 



M arcsm 



arcsm 



cos7\/l — V? 
uv — sin 7 



+ f arcsin 



u — V sm 7 



cos 



(5.6) 



Note that reduces to a in (3.14) when 7 = 0, i.e., when = Y . From (3.10), 



vr 



A(arccosM, arccosf; — — 7) = (1 — m — f )7r + 47ra^(M, t;). (5.7) 
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Figure 8. The density /■y{u,v) (with 7 = n/i) 
Next, extend the definition of a^{u, v) to C2 \ as follows (see Figure [9]): 



a^(M, v) 



' u+v—l 
4 

u—v+1 



{C2\E.f)r]{{u,v) \ u + v> l + sin7}, 
(C2 \ -E^) n {{u,v) \ v -u> I - sin7}, 
(C2 \ -E^) n {{u,v) \ v - u < sin7 - 1}, 
(C2 \ E^) n {{u,v) \ u + v < -sin7 - 1}. 

(5.8) 

Note that (5.6) and (5.8) agree on dE-^, i.e., when u^ + f^ — 2uf sin7 = cos^7. Also 



if {u,v) e i?5(7) 
if iu,v) G Rei-f) 

=^ if {u,v) e Rsii) 



for the proof of Theorem 5.3 



note that (5.8) reduces to a in (3.16) when 7 = 0. The following lemma will be useful 



Lemma 5.2. Let {U,V) be a bivariate random vector in C2 with uniform[— 1,1] 
marginals that satisfies {U,V) = {—U,—V). Then the cdf F{u,v) satisfies 

u + V 



F{u,v) 

Proof. By the symmetry condition, 



+ F{-u,-v), {u,v)eC2. 



(5.9) 



F{u,v) 



P[-U <u,-V< v] = P[U >-u,V> 

1-P[U < -u] - P[V < -v] + P[U <-u,V < -v] 

U + V 



2 J 

+ F{-u 
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□ 




Figure 9. Eight regions -Ri(7) — Rsil) ioi an elliptical copula (with 

7 = tt/S) 



Theorem 5.3. The cdf = copula of {U, Vy) is given by 

u + V + 1 



4 



+ a^{u, v), 



(u,v) G Co. 



(5.10) 



Proof. To find F^{u, v) we again use the formula (3.10) for the area of the intersection 
of two spherical caps on dB^. Here, unlike (3.12), the axes of the two caps are not 



necessarily perpendicular. The single formula (5.10) is obtained by considering the 
partition C2 = Uf^;^-Rj(7), where R5{'y) — Rs{"y) are defined in (5.8) and (see Figurejo]) 



R2{l) 

Rz{i) 
R^{i) 



n {{u,v) \ <u,v <l}, 
n {{u,v) \ -1 < u < < V < 1}, 
E^r]{{u,v) I -I <v <0 <u <1}, 
E^ n {{u,v) \ -l<u,v < 0}. 



Case 1: {u,v) G Ri{^). By using Figure 10 
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Figure 10. The region [X > ■u,Xsin(7) + y cos(7) > v] for Case 1 
(with 7 = 7r/8) 



P[U <u,V^< v] 

l-P[U>u]- P[V^ >v] + P[U >u,V^> v] 

1 - (^-y—^ - (^"~2~") ^ ^^^^^ + y COS7 > v] 

U + V 1 

1 Aiarccosu, arccost;; 7r/2 — 7) 

2 47r ^ ' ^ / !J 



U + V + 1 



+ a^{u, v) 



[by Kn] 
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Case 2: {u,v)) G R2{l)- Because {X,Y) = (— X, F) and using Figure [ll] 



F^{u, v) = P[U <u]- P[U <u,V^> v] 
u + 1 

— P[X < u, X sin 7 + y cos 7 > f ] 



2 

u + 1 
2 

u + 1 



— P[—X < u, —X sin 7 + y cos 7 > f ] 

- P[X > -u, Xsin(-7) + Fcos(-7) > v] 



2 

u + 1 1 

— — A(arccos(— m), arccosf ; 7r/2 + 7) 

2 4:71 

u + V + 1 



4 

U + V + 1 



-a_^{-u,v) [by (5.7) 



+ ay{u,v) [by (5.6)] 




Figure 11. The region [X < u,Xsin(7) + Fcos(7) > v] for Case 2 
(with 7 = 7r/8) 
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Case 3: {u,v) G Rsi'j)- Then {—u, —v) G -^2(7), so by Lemma 5.2 and Case 2, 



+ Fy{-u, -V 



u + V —u — f + 1 

H ^ h q;^(-m, -f J 



4 



U + V + 1 



+ f ) 



[by (5.6) 



Case 4: {u,v) G R^i'y)- Then (— m, — 1>) G -Ri(7), so by Lemma 5.2 and Case 1, the 
argument for Case 3 applies verbatim. 



Case 5: {u,v) G R^i'y). 



Case 6: {u,v) G Reil)- 



1-P[U >u]- P[V^ > V 
, 1 — u\ f 1 — V 

1 



U + V + 1 



P[U < u] 
u + l 



+ a^{u, v) 



[by (5.8)]. 



u + V + 1 



+ v) 



[by (5.8)]. 



Case 7: {u,v) G Rji^). Then (— u, — f ) G Rei'y), so by Lemma 5.2 and Case 6, the 
argument for Case 3 applies verbatim. 



Case 8: {u,v) G Rsil)- Then (— u, — f ) G R^i'y), so by Lemma 5.2 and Case 5, the 
argument for Case 3 applies verbatim. □ 



6. Copulas Derived from the Uniform Distribution on the Unit Ball 

Up to now we have addressed the question of whether copulas can be generated by 
means of linear functions of a circularly symmetric or spherically symmetric random 
vector. Now we ask whether non-linear functions of such random vectors can generate 
copulas. We shall restrict attention to random vectors uniformly distributed over the 
unit ball and produce relatively simple non-linear functions that generate copulas 
on Cd- 
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Figure 12. The copula F^{u,v) (with 7 = tt/A) 



We begin with the bivariate case. Suppose that (X, Y) is distributed uniformly on 
the unit disk B2 = {{x,y) e R'^ \ + < 1}. Because 



X I F ~ uniform[-Vl - Y^, Vl - Y^] 
and F I X ~ uniform[-Vl - X2, Vl -X2], 



it follows that the random variables 



satisfy 



X Y 
U:= , V:= , (6.1) 

Vl-Y^ Vl -X2 ^ ^ 



U\Y ~ uniform[— 1, 1] 



y I X ~ uniform[— 1, 1]. 
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Thus, U and Y are independent, V and X are independent, and unconditionally, 

U ~ uniform [—1, 1] 

V ~ uniform [— 1, 1], 

so the joint distribution of (f/, V) generates a copula v) on the centered cube 
C2 = [—1,1]^. Note that U and V are not linear functions of (X, Y). 

Question 6: Are U and V independent, and if not, what is the nature of their 
dependence? 

Answer 6. Clearly U and V are uncorrelated, since E{U) = E{V) = and 

XY 



E{UV) = E 



v/(l-X2)(l-F2) 



0, 



all by the circular symmetry of {X,Y). However, the joint pdf and cdf of {U,V) 
derived below show that they are not independent. 



Proposition 6.1. The joint density of {U,V) is given by 



l ^(l-^2)(l_^2) 

/K^ = . 

TT (1 — M^f "^j^ 



(6.2) 



/ 1.( 




Figure 13. Joint density f{u,v) of {U,V) in Proposition 6.1 



Proof. This pdf is again obtained via the Jacobian method. It follows from (6.1 ) that 



and 
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Substitution of the second expression for y"^ into the left side of the first relation and 
vice-versa yields 



X 



u\l-v^) 



y 



so, since x and u {y and v) have the same signs by (6.1), we obtain 



X = x[u, V) 



u 



2^,2 



y = y{u,v) 



Thus 



dx 
du 



+ u{l-u'v^)-^'\uv^) 



y2y2 



1 — u^v'^ 



dx 
dv 



u 



(1 -M2t.2)3/2' 

(l-t;2)-i/2(_^) 



2-7)2 



uv 



y2y2 



-1 H ^ 

1 — u^v"^ 



uv{l — ti^) 



By symmetry it follows that the Jacobian is given by 

Vl-t)^ uv(l-v?) 



J 



(l_„2„2)3/2 
uv{l—v'^) 



V^T^(l-«2„2)3/2 
(l-«2t,2)3/2 



V^T3^(l_u2„2)3/2 

and hence the determinant of J is given by 

^{l-u^){l-v^) 



\J\ 



(1 — u^v'^y 



Because the pdf of {X,Y) is f{x,y) = ^lB2{x,y), the result (6.2) follows. 
For < u,v < 1, {u,v) 7^ (1, 1), let Ei{u) and E2{y) be the ellipses 

{x,y) 



EAu) 



Eo(v) 



+ r <i 



[x,y) 



2 , y ^1 
x^ + — <i 



The next lemma leads to the cdf F(u,v) corresponding to the pdf (6.2). 

28 
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(0,0), 

ix{u,v), y{u,v)) 



u\/l — v'^ V\/l — V? 



2-7 )2 



2/1,2 / 



h = (n,0), 

c = (0,1;), 

d = (a/1 - y^iu,v), yiu,v)), 

e = {x{u,v), ^/l - x'^{u,v)), 

f = (1,0), 

9 = (0,1). 
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Then 



-^ATea.{Ei{u) n E2{v)) = Area(oa6) + Area(oac) 

= u Area^odf) + v Area^oeg) 

IL V 

= — arcsin(?/(u, v)) + - arcsin(a;(u, f )), 



from which the result follows. 



□ 



Theorem 6.3. The copula (= cdf) corresponding to the pdf (6.2) is given by 

, {u,v) e C2. 



F(u,v) 



U + V + 1 ^ U . ( — V? \ V . ( Ma/I — 



\ arcsm 

27r 



H arcsm 



Vl - u^v'^ ) 27r Vyr^ 




Figure 15. Nonlinear transformation copula F{u,v) in Theorem 6.3 
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Proof. Because {U,V) is sign-change invariant and has uniform[— 1, 1] marginals, it 
follows from (3.5) and (3.7) in Lemma 3.1 and from (6.1) that for {u,v) G C2, 



F(u,v) 



U + V + 1 



a{uv) P[0 < f/ < < 1/ < \v\ 



u + V + 1 aiuv] 



4 


+ 


4 


U + V + 1 


+ 


a{uv) 


4 


4 


U + V + 1 


+ 


a{uv) 


4 


Air 



P[{X,Y)eE^{u)r^E2{v)] 
Area(Ei(|M|)n ^2(1^1)). 



The result now follows from Lemma 16.21 



□ 



Problem: The construction (6.1) extends readily to generate a copula on Cd- For 
d = 3, for example, let (X, Y, Z) be uniformly distributed on the unit ball and 
define 

X _ Y Z 



U :-- 



V :-- 



W :-- 



Vl-Y^-Z^ Vl-X^-Z^ x/1 - X2 - F2 

Then the marginal distributions of U, V, and W are each uniform [—1, 1] so the cdf 
G{u, V, w) is a copula on C3. To find this copula one would need to determine 

Volume (El (u) n E2{v) n E^iw)), 

where now, for < u,v,w < 1, Ei{u), E2{v), and E^{w) are the ellipsoids 



Ei{u) 
E2{v) 
Esiw) 



{x,y,z) 
{x,y,z) 



X 



+ y <i 



2 . y I 

X + — + z < 1 



„2 I „.2 _|_ 



+ + ^ <l 
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